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Abstract 

We present a model for the structure of baryons in which the valence partons 
interact through a linear potential. This model can be derived from QCD in the 
approximation where transverse momenta are ignored. We compare the baryon 
structure function predicted by our model with a standard global fit to Deep 
Inelastic Scattering data. The only parameter we can adjust is the fraction of 
the baryon momentum carried by the valence partons. Our prediction agrees 
well with data except for small values of the Bjorken scaling variable. 
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The quark-parton model was proposed as a simple explanation of the 
phenomenon of scaling in Deep Inelastic Scattering. After the discovery of 
Quantum Chromodynamics (QCD), scaling was understood as a consequence 
of asymptotic freedom Moreover, QCD predicts small violations of scaling 
that can be calculated perturbatively: knowing the parton distribution at one 
(large enough) energy scale, it is possible to predict them at any other large 
energy. Although these predictions have been confirmed impressively, there 
is as yet no theoretical understanding of the initial parton distributions; they 
are controlled by non-perturbative effects. Accurate knowledge of the parton 
distributions is essential in order to make predictions about any high energy 
scattering process in a hadron collider. Also, only after we understand the par- 
ton distributions as predicted by the standard model, can we look for deviations 
which may signal a substructure for quarks. As a result tremendous effort has 
been expended by several groups of theorists and experimentalists to extract 
the parton distributions from scattering data at various energy scales ^, |^] . 

We will present a model of interacting partons for hadronic structure func- 
tions. The partons are assumed to be relativistic particles interacting with 
each other through a linear potential, which is a sensible idea in the lightcone 
formalism. Their momenta transverse to the direction of the collision will be ig- 
nored. We will solve for the parton wavefunction within the approximation that 
it factorizes into a product of single particle wavefunctions. This is analogous 
to the Hartree-Fock approximation of atomic physics. After these approxima- 
tions, the wavefunction can be determined as the solution of a nonlinear integral 
equation. This equation can be solved numerically, although we also find a good 
approximate solution by analytic methods in some special cases. 

We can compare with the structure functions extracted by several groups of 
experimentalists and theorists. Considering that we have at most one parameter 
to adjust, the agreement is quite good, except for small values of the Bjorken 
scaling variable xb- In this region our model is not expected to be valid: sea 
quarks and gluons cannot be ignored. 

Can we derive this interacting parton model from QCD? In fact we have al- 
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ready given such a derivation in previous publications Q . In the approximation 
in which transverse momenta are ignored, two of the space-time co-ordinates 
become irrelevant: QCD is replaced by its dimensional reduction to two dimen- 
sions. That such a 'coUinear QCD' can describe hadronic structure functions 
has also been proposed by other scientists Q] . What is unique to our approach is 
the picture of a baryon as a topological soliton, which remarkably enough, gives 
a derivation of our interacting parton model from first principles. In previous 
work, two dimensional QCD was shown in the large N limit to be equivalent 
to a bilocal field theory whose phase space is an infinite dimensional Grass- 
mannian. The baryon is a topological soliton in this theory, whose energy and 
structure can be estimated (within a variational approximation) by a nonlinear 
integral equation. This is precisely the integral equation we will study here. 
After the equation was derived and studied in this way it became clear that 
it had a simple interpretation in terms of the parton model. It is this parton 
model point of view that we will mostly use in this paper. 

In this approximation, the dynamical degrees of freedom associated with 
gluons are ignored: only the longitudinal components of the gluons which can 
be eliminated in favor of the quarks are kept. In later publications we will study 
how tranvcrse gluons will modify the structure of a baryon; in particular derive 
the gluonic distribution functions. 

The stucture functions of a hadron are measured in deep inelastic collisions 
of a virtual photon ( or weak gauge boson) with a proton or neutron. Since 
the component of momentum in the direction of collision is much greater than 
those transverse, it is usually assumed that all the constituents of the hadron 
(partons) are moving in the same direction. We will choose this direction to 
be the x^-axis. It is also customary to use the null component of momentum 
p ~ Po — Pi instead of the spatial component as the basic kinematic variable 

H . In terms of the energy Pq and null momentum, the mass shell condition 
for a particle of mass m is 

- p^ + 2pop = 



3 



or equivalently, 

Note that since po — ^{p\ + m?) > pi, the null component of momentum 
p ~ Po ~ Pi is always positive. ^ 

Let ^^{v^ ct,p) be the wavefunction of a single parton expressed as a function 
of the null component of momentum. Here a = 1, • • • M is a discrete quantum 
number that labels spin and flavor; and u = 1, ■ ■ ■ N labels color. Although the 
value of color in nature is three, we will find it convenient to keep it arbitrary 
for now. Because of the positivity condition on p, 

ip{i^, a,p) = 0, for p < 0. 

A baryon, which is made up of N such partons, will have a wavefunction 
•0(^1, ai,pi; 1^2, 012, P2', • • • ; i^N , cuntPn)- This will vanish unless all the momenta 
are positive. Since the partons are fermions (after all, they are quarks), the 
wavefunction must be totally anti-symmetric under permutations of the par- 
tons. However, the baryon must be colorless so the color quantum numbers 
must be completely antisymmetric by themselves: 

ii{vi,ai,pi;v2, a2,P2; - ■■ , vn^un^Pn) = <^vt.M2-r--'^N'^{oii,pi]a2,P2] ■ ■ ■ cun^Pn) 

where Cui,v2,---un Levi-Civita tensor. So we can forget about color and 

deal with the remaining part of the wavefunction ■0(q;i,pi; a2,P2; • • • ; olntPn) 
which is completely symmetric. 

The kinetic energy of a system of N partons will be, 

^' 1 2 J 7 

Ir , "^a, II 7/ x,2dpi ■ --dpN 



N 



•^The mass shell condition looks a bit unusual bcacusc the metric tensor is not diagonal in 
our choice of variables. If we use the conjugate variables (it, x) to po,P, ds^ = du^ + 2dudx. 
Here x is the usual spatial co-ordinate and m = — a; is a null co-ordinate. 
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If there is a two body potential g^w(x)between the partons,the total potential 
energy will be Q 



1 r°° 



Xi — Xj)\ip{ai,xi\ • • • ; ajv, XN)\^dxi ■ ■ ■ dx 



N 



Here, 

ip{ai,xi;- ■ ■ ai,XN) ^ I tp[ai,pi;- ■ ■ aN,PN)e — 7o~Tn — 
Jo 

is the wavefunction in position space. 

In addition to the kinetic and potential energies, there could be a term 
in the hamiltonian describing the self-energy of the partons. Due to Lorentz 
invariance, such a term can only be a (finite) renormalization of m^; i.e., it is 
of the form 

N 



roc ^' g2 

E / E^i^("i'Pi;---;"^' 



..2dpi ■ --dpN 



for some constant c. 

Thus, the ground state wavefunction is determined by minimizing the total 
energy 

^ 1 ,,2 



2 dpi • --dpN 
Pi; - ■ ■ aN,PN)\ — tkzTn — 



(2^) 

+ 2^^ E / y^^(a:» - a:j)|'0(ai,xi; ■ ■ ■ aN,XN)\'^dxi ■ ■ -dxN- 

Here /i^ = 77?? + is the renormalized mass of the parton. 

The success of the simple parton model and its QCD inspired variants shows 
that it is a good approximation to describe the system in terms of particle distri- 
butions that depend only on the momentum of a single parton. In other words 
correlations between the partons can be ignored, except in that the total mo- 
mentum of the partons must be fixed. Thus it should be a good approximation 

A multiplicative factor has been chosen such that the coupling constant g has units of 
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to assume that (in the ground state), the baryon wavefunction has the form : 

N 

'(/'(ai,Pi; • • ■aN,PN) = 2ttS{P - y^Pi) ]~[ -0(ai,Pi). 

i i—1 

Here, P is the momentum of the baryon. Since each of the parton momenta pi 
are positive, we see that 

ip{p) — unless < p < P. 

This ansatz is anafogous to the Hartree approximation of atomic physics, 
which works weh even when the number of electrons is small as in the Helium 
atom. In any case we can regard our product as a variational ansatz for the 
ground state of the baryon. We will normalize the single parton wavefunction 
to have length one: 

Moreover it must satisfy the sum rule on momentum: the total momentum of 
all the partons must equal the momentum of the baryon: 

N j\\4^{p)\'^=P 

In fact the valence partons may not carry all the momentum of the baryon. 
(The structure functions extracted from data show that only about half the 
momentum of the baryon is carried by the valence partons.) We can allow for 
this by replacing the momentum sum rule above by 

N r p\m?^^fp 

Jo 

where / is the fraction of the baryon momentum carried by all the valence 
partons. We will see that this fraction / is the only parameter on which the 
sturcture functions will depend. 
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The energy per parton becomes now 

1 1'°° 

(1) -~g' z;(x-2/)^|V'(a,x)p^|V'(/3,y)pdxdy. 

where, g'^ ^ g'^N and ^{a,x) = ip{p)e'-P'' ^. 

What potential v{x) should we use? It is known that the potential between 
quarks is, to a good approximation, linear. Moreover the coUinear approxima- 
tion to QCD mentioned above predicts a linear potential. Hence we will choose 

A linearly increasing potential can lead to an infrared divergence in the energy of 
a wavefunction, even one that is decaying at infinity. For example if "0(0!, x) ~ ^ 
for large x the potential energy term will diverge logarithmically. This is where 
a proper choice of the self-energy term mentioned above comes in. If we choose 
the constant c in the self-energy such that 

2 2 ff^ 

Ma = "^a : 

TT 

this Infrared divergence dissappears. This is indeed the value of the self-energy 
term predicted by ( the large iV limit of ) two dimensional QCD. 

The wavefunction of the partons in a baryon is now determined by minimiz- 
ing the above energy function subject to the constraints on the wavefunction 

0(p) = for p < andfor p> P and | |^ — 1- 

The minimization of the above energy is equivalent to the solution of a 
system of nonlinear integral equations. The condition that the wavefunction 
vanish outside of the interval [0, P] in momentum space is a subtle constraint 
on the position space wavefunction: 7/1(0,0;) must be the boundary value of an 
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entire function which grows at most hke e^'^' at infinity. Since such a condition 
is clearly difficult to impose numerically, we work in momentum space. Variation 
of the energy with respect to ip gives 

1 1 d 



where 



^ dq 
2tt' 



P ^ Jo 



Here, A is the Lagrange multiplier that enforces the constraint = 1- (Since 

the function we are minimizing is quartic in this is not the same as the 
minimum energy E.) 

The second equation is just the momentum space version of the equation for 
the mean potential due to all the partons: 



V"{x) = \^{x)\ 
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Note that the Fourier transform V{p) = V{x)e ^"P'^dx of the mean potential 
is singular at the origin: 

Hence the integrand in the equation for ^p{a,p) is singular. We must define it 
to be a finite part integral in the sense of Hadamard : 

V{p-q)ipia,q)— ^ / [ijj{a,p + q) + 4j{a,p - q) - 2ip{a,p)]V{q) — . 



By studying the behavior of the equation in the neighborhood of the singular 
point p ^ 0, with the ansatz ip{p) ^ p'^ , we get the following formula relating 
the exponent v and m^: 

Thus, at the critical point to = 0, we have = 0: the critical wavefunction 
tends to a constant as p ^ 0. If m? > 0, the wavefunction vanishes like a 
power: v > 0. 



9' 



1+ ^ 



dy 

2 



y 



{i+yr + {i-yr 



dy 

,,2 
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This nonlinear integral equation was derived previously ||^ from the large 
iV-limit of two dimensional QCD. The 'master field' of this limit of 2DQCD is 
a hermitean matrix M{a,p; [3, q) satisfying the nonlinear constraint 

/dr 
— M{a,p;-f,r)M{-f,r;/3,q)^0. 

The set of solutions of this condition is an infinite dimensional Grassmannian 
manifold, which has connected components labelled by an integer: the 'vir- 
tual rank', —5X^0/ ^^{(^tP\'^iP)^- This integer has the physical meaning of 
baryon number. The baryon is thus a topological soliton and is described by a 
configuration of minimal energy in the sector of virtual rank one. 

Thus we have here a realization of Skyrme's idea that the baryon is a topo- 
logical soliton 1^. The energy of the configuration M{p, q) is given by 
,2 



m) = Yl fl{p+—)M{a,p;a,p)^+9'Y. f \M{a,x;p,y)\ 



'v{x — y)dxdy. 



Here, 



M{a,x;P,y) = / M(a,p; /3, g)e'f--*«^^. 



(2^) 

The parameters p^ and g are related to those of 2DQCD: 

where m is the current quark mass and g the gauge coupling constant. Note that 
in the limit of chiral symmetry, rua = 0, (when the lightest meson is massless) 
and the value of /i^ is actually negative: 

This is the value we will mostly study. 

In addition to topological soliton solutions describing the baryon the above 
equation also describes small fiuctuations from the vacuum. The corresponding 
particles are the mesons. If we assume that Af(p, q) is infinitesimally small, its 
equations of motion will reduce to a linear integral equation. This equation 
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was first derived by 't Hooft by a masterful use of diagrammatic methods pO| . 
Witten suggested later that the baryon can be described by a Hartree-Fock 
aproximation in the large N limit of QCD. He carried out this idea in a non- 
relativistic context. Our model may be thought of a relativistic implementation 
of this idea of Witten. 

Now if we make the variational ansatz that M{a,p; f3, q) = —2i/j{a,p)i/j* {(3, q) 
is separable, the quadratic constraint on M becomes the condition that 

ip(a,p) — ior p < 0. 

Moreover, the energy of the soliton reduces to just the formula we obtained 
within the parton model! (The momentum sum rule is the the analogue of the 
WKB quantization condition in this theory.) This is how we can derive the 
parton model from coUincar QCD. It is remarkable that the topological soliton 
model, which at first looks so different from the parton model, leads to exactly 
the same equation for the hadron structure function. 

The variational approximation that M be separable corresponds in the par- 
ton model to ignoring the 'sea quarks'. This seems reasonable only for values 
of the Bjorken parameter xb that are not too small. In general we should get 
better answers by directly minimizing the energy of the soliton model without 
the further approximation of the separable ansatz. We are currently studying 
this issue. 

Let us return to solving the integral equation. It cannot be solved ana- 
lytically, so we will resort to a numerical method. ( We can get approximate 
analytical solutions in some special cases which agree well with the numerical 
solution.) For simplicity, we will from now on put all the parton masses equal to 
each other. It is not difficult to modify our solutions to include unequal masses. 
Also, it is sufficient to find a solution that is non-zero only for one value of the 
spin-flavor index a: 

■4'{a,p) = (5a,iV'(p). 

This breaks the U{M) invariance of our model spontaneously. This symme- 
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try can be restored later by the collective variable method as in the theory of 
solitons, but we will not address this issue here. 

Let us first describe our approximate analytic solution for the special case 

~2 

^ = and large N. We are mostly interested in the case /i = — ^ but 
this approximate solution will help to validate our numerical procedure. It is 
reasonable to expect that most of the contribution to the integral 

mp + q) + ^(p - q) - 2i;{p)]V{q)^ 

Ztt 

will come from the neighborhood of the singularity near q — Q: 

\i'{p + q)+ ^{P - 9) - 2i;{p)]V{q)^ ^ r{p) r eV{q)^. 

^TT Jo 



Set 



r 2t>, 
Jo 



Using the fact that the ground state wavefnction is real, we can show that 



Thus we get a differential equation: 

-fai^"ip) + l[p+^-2X]i^ip)^0. 
2 P 

This can be thought of as a Schrodinger equation with a linear potential plus a 
Coulomb potential; such an equation has been used to study heavy quark-anti- 
quark bound states [^2j. But the physical origin of our equation is completely 
different: we get this equation in momentum space not position space. Also, 
our quarks are not assumed to be heavy and we are studying a baryon not a 
meson. 

In the special case /j,^ = 0, the solution is a linear combination of Airy 
functions. But it is awkward to work with this linear combination. If we make 
the further approximation that N is large, our momentum sum rule shows that 
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P is large as well. Then the boundary condition on the wavefunction is that it 
must tend to zero at infinity: in other words 



The eigenvalue A is determined by the condition that the solution vanish at 
the origin: 

Here, ~ —2.33811 is the root of the Airy function closest to the origin. C is 
determined by the normalization condition. 

The constant a is determined by putting this solution back into the definition 
of a. We get, 

2a = g 3- 3-, A ^ „ 7, ' , ^ ~ 0.847589g. 

(2-)^ [/- Ai ^ ^(2-) [J- Ai ^{Odi] ^ 



Moreover 



^(p) = CAi (ei(^-l))- 



This analytic approximation agrees well with the numerical solution (de- 
scribed below) for /i^ = and large N . This confirms the validity of our nu- 
merical method. 

Now let us discuss the numerical solution of this problem. This is not 
straightforward since the kernel of the integral equation is singular. We need a 
reliable method of numerical quadrature for integrals such as 

V I f(.P,q)p{q)dq 
Jo 

when the weight function p{q) has a singularity at g = like p-. We need 
to subdivide the interval [0, P] = U"^]^[6r, &r+i] into subintervals. Within each 
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subinterval we choose a set of points qjr,j = l,'''^^r- We approximate the 
integral by a sum 

rP 

/ fip, q)p{q)dq = ^ Wjrf(qjr)- 



The weights Wjr are determined by the condition that within each subinterval 
\br, fer+i], the integral of a polynomial of order — 1 is reproduced exactly: 

V I q^p{q)dq = 2_, ^jrq^. 

This is just the usual method of numerical quadrature except that the integral 
on the left hand side is singular for hr — Q and fc = 0, 1. In these cases we 
can evaluate the l.h.s analytically as the finite part in the sense of Hadamard. 
(The main difference from the usual situation is that the moments on the l.h.s 
of the above equation are not all positive.) The weights are then determined by 
solving the above system of linear equations. 

Given an approximate mean potential Vs{p) we can convert the linear inte- 
gral equation 

-[p+^- 2A,]^,+i(p) +3^7' / V,{p-q)A+i{q)Tr = « 
2 p Jo 2-K 

into a matrix eigenvalue problem by the above method of quadrature. We use 
the ground state eigenfunction so determined to calculate numerically the next 
approximation Vs+\ for the mean potential. This process is iterated until the 
solution converges. Having determined the wavefunction, we must impose the 
momentum sum rule to determine g. 

We chose n = 50 intervals each containing 1/^ = 6 points for our quadrature. 
The resulting matrix was diagonalized using Mathematica. As a starting mean 
potential we used 




The procedure converges in about 10 iterations. We find that due to the errors 
in the discretization of the problem, the critical value of /^^ (the value at which 



13 



~2 

the wavefunction goes to a constant at p = 0) is a bit lower (— 1.3pj) than the 

~2 

theoretical value of — — . 

Now we turn to the question of the comparison of our model with data from 
Deep Inelastic Scattering. It is customary to describe the parton distributions 
as a function of the Bjorken scaling variable < < 1 which is the fraction 
of the null component of momentum carried by each parton. This means we 
will rescale momenta to the dimensionless variable p = xbP- The probability 
density of a parton carrying a fraction xb of the momentum is then 

(The factor of is needed because (j){xB) is traditionally normalized to one 
with the measure dxB rather than 

It is important to note that the only dimensional parameter in our theory, g, 
cancels out of the formula for (^{xb)'- it only serves to set the scale of momentum 
and when the wavefunction is expressed in terms of the dimensionless variable 

it cancels out. We have set /i^ to the critical value (within nmiierical errors), 
which is the value corresponding to chiral symmetry; i.e., zero current quark 
mass. The number of colors we fix at = 3. Thus the only free parameter 
in our theory is the fraction / of the baryon momentum carried by the valence 
partons. The parameters N and / appear in the combination iVoff = y. 

We have ignored the ispospin of the quarks in the above discussion. We 
should therefore compare our structure functions with the isoscalar combination 
of the valence quark distributions of a baryon, 0(xb). It is not difficult to take 
into account of isospin effects. 

The parton distributions have been extracted from scattering data by sev- 
eral groups of physicists In the accompanying figures we plot our 
wavefunctions for a few values of / to those extracted from data by the MRST 
collaboration. We agree remarkably well with experiment except for small val- 
ues oi Xb- The agreement is best when the fraction of the baryon momentum 
carried by the valence partons is about / — 0.6. 

Our model does not predict the observed behavior of the parton distributions 



14 



for small xb' our probability distribution tends to a constant for small xb 
although due to numerical errors that is not evident in the plots. The observed 
distributions have an integrable singularity there: roughly speaking, 4>{xb) ~ 
x^'^ for small x^.The approximations we made clearly break down in the small 
Xb region: sea quarks can no longer be ignored, indeed even gluons need to be 
considered. We will study these effects in future publications. 
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Numerical and Analytic Wave Functions 
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Comparison of parton wavefunctions \/(j>{x) with the MRST global fit to data. 
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